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We calculate certain estimates for the solution of the characteristic problem of the wave equation 
reduced to first order, in terms of the free data prescribed on two transverse surfaces, one of which 
is characteristic. Estimates of such kind ensure the stability of the solutions under small variations 
of the data. Similar estimates exist for the derivatives of the solution as well. 
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WELL POSEDNESS AND CHARACTERISTIC 
PROBLEMS 



Given a system of partial differential equations where 
a unique solution exists for some given data, oftentimes it 
is necessary to know how sensitive the solution is to small 
variations of the data. One would not want the equations 
to amplify any uncertainty in the data beyond control. 
When, in addition to existence and uniqueness, the equa- 
tions guarantee the stability of the solution under small 
perturbations of the data, the problem is referred to as 
well posed p|. 

There are two components to a well-posed problem. 
First there are the equations themselves, which are usu- 
ally defined only up to transformations of variables and 
coordinates - respecting the differential order of the equa- 
tions. Second, there is the data set, usually an appro- 
priate set of values prescribed on a chosen surface. In 
coordinate-independent terms, the choice of data surface 
determines the type of problem for the given equations. 
If one seeks a solution in the interior of a region bounded 
by a closed data surface, we have a boundary problem. If 
the solution is needed on half of the available space of 
independent variables (x 1 > in the appropriate choice 
of coordinates for the problem), then we have an initial 
value problem, or "Cauchy problem" . It should come as 
no surprise that some equations have at the same time a 
well-posed problem and an ill-posed problem. The text- 
book example is the two-dimensional Laplace equation 
u xx + u yy = 0, for which the boundary problem is well 
posed but the initial value problem is not. A less acknowl- 
edged but equally powerful example is the wave equation 
Uxx — Utt — 0, for which the reverse statement holds. 
Typically elliptic equations admit well-posed boundary 
problems, whereas hyperbolic equations admit well-posed 
initial value problems. 

In the case of hyperbolic equations, there are also char- 
acteristic problems. In this case, for a unique solution to 
exist, the data must be prescribed on two intersecting 
transverse surfaces, one of which is ruled by character- 
istic curves and we assume the other one is not Q- On 
the side, characteristic problems are not to be confused 
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with mixed initial-boundary value problems, where val- 
ues are prescribed on two intersecting surfaces but such 
that the initial surface is not characteristic, because in 
such mixed problems the number of data variables nec- 
essary to specify a unique solution is larger than in the 
initial value problem of the same equations. The charac- 
teristic problem of a set of equations requires the same 
number of data as the corresponding initial value prob- 
lem. In characteristic problems, the method of solution 
is an adaptation by Duff of the Cauchy-Kowalewsky the- 
orem for initial value problems. As for the stability of the 
solutions, there does not seem to have been a great deal 
of interest. Characteristic problems are usually neglected 
in favor of their very powerful cousins, the initial value 
problems, which require only one correctly chosen data 
surface. Still, in some physical contexts such as general 
relativity, characteristic problems have traditionally been 
very fruitful in connection with issues of radiation 
and have been used methodically to obtain numerical so- 
lutions Q . The question of well-posedness in such cases 
becomes quite relevant. 

Here the stability of the solutions of the characteristic 
problem of the wave equation in three spatial dimensions, 
reduced to first order, is explored by means of a partic- 
ular method. The method is analogous to the energy 
estimates used in the case of initial value problems of 
first-order hyperbolic systems 0, 0. The principle be- 
hind the method is to obtain, as a consequence of the 
equations, an inequality by which the "size" of the so- 
lution is bounded by the "size" of the data, namely: to 
estimate the solution in terms of the free data. In linear 
problems, the estimate applied to the difference between 
two solutions that are initially close trivially ensures sta- 
bility. 

In Section[n]we include a brief expository review of the 
method for the initial value problem of the wave equation 
reduced to first order. This is done with the purpose of 
developing a parallel with either characteristic planes or 
cones as data surfaces in Sections IIIII and llVl It is found 
that a certain kind of estimate for the solution in terms 
of the data can be established in both cases, demonstrat- 
ing that these two characteristic problems for the wave 
equation are well posed. Section Ivl contains concluding 
remarks including a brief overview of the status of the 
question of stability of characteristic problems in a more 
general sense. We hope to uphold that even classic prob- 
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lems such as the wave equation have the potential for 
new insights when viewed from a different angle. 



II. THE CAUCHY PROBLEM OF THE WAVE 
EQUATION REDUCED TO FIRST ORDER 

We have the wave equation in cartesian coordinates in 
three spatial dimensions given by 



i>tt = 4>xx + Ipyy + Ipz 



(1) 



where a subscript x a = (t, x, y, z) denotes partial differ- 
entiation with respect to the coordinate usual. We 
do not deal with this equation directly, but prefer to "re- 
duce" (actually, extend) the equation to a system of first 
order equations. In doing so, one substitutes the original 
problem with a problem that has more solutions, among 
which the solutions to the original problem can be singled 
out. 

We define new variables U = ipt > P = ipx > Q = ipy an d 
R = ipg. In terms of these variables, is equivalent to 
the following system 



Pi 

Qi 
Rt 

It 



U x , 

Uy, 
U 2 . 

u, 



(2a) 
(2b) 
(2c) 
(2d) 
(2e) 



in the sense that all the solutions ij){x a ) of JIJ are so- 
lutions of J3J. In order to single them out we need to 
impose the following constraints on the initial data sur- 
face 



C X =P-^ X = 0, 
C v = Q-^ y = 0, 
C z = R-ip z = 0, 



(3a) 
(3b) 
(3c) 



which remain satisfied because their t— derivatives vanish 
(C x t = C y t = C z t = 0) by virtue of ©. However, we do 
not have to restrict the data for the estimates that follow, 
so we ignore the constraints from now on. In matrix 
notation the system @ has the form 



A a d a v + Dv = 



(4) 



where v — (U, P, Q, R, ip), A* is the identity matrix, the 
other three matrices A 1 are symmetric, D has all van- 
ishing coefficients but one (of value —1), and summation 
over repeated indices is understood. Multiplying on the 
left with v we have 



vA a d a v + vDv = 0, 



(5) 



and, because the matrices are symmetric, we can extract 
the partial derivatives to obtain a "conservation law" : 



Because we are not interested in a boundary-value prob- 
lem, we now assume that the fields either decay suffi- 
ciently fast at infinity so that their volume integration in 
the entire space R 3 is finite, or else they are periodic, so 
that we can restrict attention to a box in R 3 . Equiva- 
lently, we assume that the fields have a Fourier transform 
or series. With this assumption we can first split off the 
time part of the conservation law © , and then integrate 
it in the spatial volume V (denoting either R 3 or a box): 



— / vA f v (Fx + / d l {vA l v) d A x + 2 / vDv d A x = 0. 

With our assumptions, the second term evaluates to zero. 
In the first term we have the Li norm of the fields, be- 
cause A* is the identity matrix and thus 



vA v d x = \ \v(t, ■ 



(8) 



Since vDv — — Ui/j, Eq. Q is equivalent to 

j t \\v(t,-)\\ 2 = J 2Ui,d 3 x. (9) 

Finally, since 2Utp < U 2 + tJj 2 then J y 2Uip d 3 x < 
\\v(t, -)||, so we have 



dt 



IKV)H 2 <IKV 



which implies 



IKMII^e'Ho, 



(10) 

(11) 



and we have our estimate of the "size" of the solution 
in terms of the "size" of the initial data. Since the sys- 
tem is linear, the difference of two solutions is also a 
solution, so small variations of the initial data result in 
small variations of the solution, at least for small enough 
times. The exponential factor is there because of the 
presence of undifferentiated terms in In spite of the 
exponential factor, the inequality shows that the solution 
depends continuously on the data, so that the variation 
in the solution at any time can be controlled by refining 
the accuracy of the data. Nonetheless, as usual the 
inequality is useless for large values of t for numerical 
purposes, where the interest resides in explicitly estimat- 
ing the error in a given solution that starts with no less 
than round-off error. 

Alternatively, one can get a better estimate in this case 
by treating ip separately from the rest of the variables. 
Clearly the equations for the variables U, P, Q and R de- 
couple from ip, and have no undifferentiated terms. A 
solution can be found independently of the value of ip, 
and can be estimated also independently of ip. If we re- 
produce the preceding reasoning with v = (U, P,Q, R), 
then we have Eq. (0} with D = 0, which leads to 



d a {vA a v) + 2vDv = Q. 



(6) 



dt 



KMII 2 = o, 



(12) 
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and consequently 

||«(V)|| 2 = IK(V)|| 2 , (13) 

without an exponential factor. Then ip can be estimated 
by 

^ = ^|t=o+ / Udt' (14) 
Jo 

which is the solution of ipt = U for known source U and 
given initial data at z = 0. This leads, for instance, to 
an estimate of the form 

ijj 2 < 2V|f =0 + 2 (J Udt'^j . (15) 

Most remarkably, one can also control the smoothness 
of the solution. The spatial derivatives of the fields have 
estimates in terms of the spatial derivatives of the ini- 
tial data. One can see this easily by taking a spatial 
derivative d{ of the system I|2al2d|) . which thus becomes 
an evolution system for the variables vi — (Ui, Pi, Qi, RP) 
of exactly the the same kind as that for v, leading to esti- 
mates for each spatial derivative Vi in terms of the same 
spatial derivative of the initial data: 

|k(V)|| 2 = lk(0,-)H 2 , (16) 

for i — x,y, z alternatively. 

As a hyperbolic equation, the wave equation admits 
characteristic lines at 45° with the time axis in all spa- 
tial directions from any fixed point, that is: all light like 
straight lines through any point. A (hyper)surface that 
has a null (lightlike) normal at all points is automati- 
cally ruled by characteristic lines and is referred to as a 
null surface. A null surface is simply the evolution of a 
wavefront forward in time out of any given initial shape. 
Any null surface can be chosen as the data surface of 
a characteristic problem. Thus there are infinite ways 
to set up a characteristic problem for the wave equation. 
But there are two aesthetically appealing cases. The first 
case is that of a null plane acting as the data surface (the 
evolution of a plane wavefront). Of course, there are null 
planes in all possible directions, so this may not be useful 
unless our problem has a particular preferred direction. If 
there is no particular preferred direction then null cones 
are suitable (the evolution of a spherical wavefront out of 
some point of origin). These two choices of characteris- 
tic data surfaces lead to similar characteristic problems, 
which are examined in the following two Sections. 

III. WELL-POSEDNESS OF THE 
CHARACTERISTIC PROBLEM OF THE WAVE 
EQUATION WITH A PLANE 
CHARACTERISTIC DATA SURFACE 

The simplest characteristic problem of the wave equa- 
tion prescribes data on a plane characteristic surface. In 



addition, we use a plane timelike surface to prescribe the 
necessary complementary data. In Subsection 1111 Al we 
examine the stability of the solutions under small per- 
turbations of the data. Estimates for the derivatives of 
the solutions are derived in Subsection IIII Bl 

A. Estimate of the solution 

We transform coordinates (t, x, y, z) — » (u, x, y, z) with 

u = t — z, (17) 

so that the level surfaces of u are null planes. The wave 
equation Q in these coordinates reads 

Zlpuz - Ipxx - 4>yy ~ ^zz = 0. (18) 

The second derivative with respect to u does not appear 
in the equation, so we only need to define three additional 
variables to reduce the equation to a first-order system. 
We have, as before, P = ip x ,Q = ip y and R = ip z . In 
terms of these variables, Q is equivalent to the following 
system 

2R U = P x +Qy+Rz, (19a) 
Pz = R x , (19b) 

Qz = Ry, (19C) 

= R, (19d) 

in the sense that all the solutions ip are contained in the 
set of solutions of this system. This system has a unique 
solution in the "wedge" space with u > and z > if 
the value of R is prescribed on the null surface u = and 
the values of P, Q and ip are prescribed on the surface 
z = 0. This is referred to as the canonical form of the 
characteristic problem of the wave equation Q . Like all 
characteristic problems in canonical form, it consists of 
two sets of equations: one of them internal to each char- 
acteristic surface [Eqs. I|19bll9djl ]. and the other moving 
from one characteristic surface to the next [Eq. I|19a|) ]. 
The field R is a normal variable, whereas P, Q and ip 
are referred to as null variables Q{. We can single out 
solutions to the original second-order wave equation by 
choosing data so that the constraints l|3a(l and l|3b|l are 
satisfied along the surface z = 0. They remain satisfied 
everywhere because d z C x = d z C v = by virtue of ifT^I) . 

As in the case of the initial value problem in the pre- 
ceding Section, we see that the equations for (R, P, Q) 
decouple from ip. Therefore we can find and estimate a 
solution independently of ip. Disregarding Eq. I|19d|) . the 
rest of the system (|19|) has the form 

A u v, u +A z v, z +A x v, x +Av v , y = (20) 

where v = (R, P, Q), and the matrices are given by 

/200\ / -1 0\ 

A"= A z = 010 (21) 

\oooy \ o o 1/ 
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A y = 


■ 





• 



















A x = 



Because all the matrices A a (a 
ric, it follows that 



(22) 



u, x, y, z) are symmet- 



O a (vA a v) = 0, 



(23) 



so we obtain a "conservation law" of the form d a J a = 0. 
Integration of this conservation law on any given volume 
in R 4 (of measure d 4 x = dudzdxdy) gives us relation- 
ships between the flows of vA a v across the surfaces that 
enclose such volume. We choose our four-volume of in- 
tegration to be a "hyper-prism" , bounded by the sur- 
faces it = (of outward pointing normal n\ = —5%), 
2 = (of outward pointing normal n 2 — —5*) and 
u + z = T for a fixed constant T (of outward pointing 
normal n\ = (5™ + 6D/V2). Assuming that the fields 
admit Fourier transforms in the x and y directions, we 
let the "hyper-prism" extend to infinity along x and y, 
or alternatively we can have a periodic boundaries in x 
and y. Our four-volume of integration V 4 is represented 
in Fig. ^ Integrating on this "hyper-prism" we find 



0= / d a (vA a v)d'x = 
/y 4 



v(A u 

u+z=T 



-A z )v dzdxdy 



vA u v dzdxdy 
vA z v dudxdy (24) 

z=0 

Since, by J2U, A u +A z is the identity matrix, then v(A u + 
A z )v = R 2 + P 2 + Q 2 and its integral on the surface T, T 
may define a (positive definite) norm of the solution, 



|2 _ 

It = 



v(A u + A z )v dzdxdy, (25) 

u+z=T 



which may be used as a measure of the "size" of the 
solution generated by data on u = and on z = 0. On 
the other hand, vA u v = 2R 2 and vA z v = -R 2 +P 2 +Q 2 . 
With these expressions, (|24f) can be written as 



Mil 



2 /i? 2 dE u + [{-R 2 + P 2 + Q 2 )dT, z 
is„ is. 

which trivially implies the following inequality 



Mlr<2 



R 2 dT, l 



P 2 + Q 2 ) tffij 



(26) 



(27) 



The right-hand side is a positive-definite measure of the 
"size" of the free data, properly taking into account both 
data surfaces. This differs markedly from other possible 
energy-like norms that could be defined, mirroring the 
Cauchy problem, by integrating the squares of all fields 
on the initial characteristic surface. When all fields are 




FIG. 1: The "hyper-prism" of integration for the characteris- 
tic problem with a null plane data surface. The top bound- 
ary surface Et lies at u + z = T constant and has measure 
dT>T = dxdydz/^/2. The surface E u is a null plane at u — 
and has measure dE u = dxdydz. The surface E z lies at z = 
and has measure dE z = dudxdy. Both surfaces E„ and T, z 
are data surfaces for the characteristic problem. 



integrated on the surface at a fixed value of u, one is, in 
a sense, over-counting the null data, thus a literal trans- 
lation of the energy of the Cauchy problem does not, in 
fact, play the intended role. On the contrary, the inequal- 
ity lf2*7|l functions as the proper characteristic analog of 
the standard estimate for initial value problems. Small 
data R on u = and P, Q on z — result in small values 
of R, P and Q within the "wedge" space u > 0, z > 0. 

Once R is known, we can estimate ip by integrating 
Eq. i|19d[) with given values on z = 



i> = ^| 2=0 + / Rdz', 



which implies, for instance, 



ip 2 < 2^|Lo + 2 ( J R dz ' 



(28) 



(29) 



B. Estimates of the derivatives 

It is not as simple to estimate the derivatives of the 
variables in terms of the derivatives of the data as it is 
in the case of the Cauchy problem discussed in the pre- 
vious Section. However, with some ingenuity, we can 
show that the derivatives R x , R v , R Zl P x , P y , P u , Q x , Q y 
and Q u can be controlled with the choice of data. These 
play a role that is analogous to that of the space deriva- 
tives in the case of the Cauchy problem discussed in Sec- 
tion[n]- P Z ,Q Z and R u are already given by the system of 
equations (|I9|) , and can be estimated directly by the use 
of the system of equations, as in the case of the Cauchy 
problem - . Of these, R x ,R y ,R z should be treated as 
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normal variables and the rest as null variables. If so, 
the data for R x , R y , R z can be obtained from the deriva- 
tives of the data for R and we avoid the introduction of 
more arbitrary data (likewise with the rest of the first- 
order variables). With this in mind, we obtain equations 
for R x ,R y ,R z by applying di to Ijl9a|l . Every appear- 
ance of P z or Q z in the result is substituted in terms of 
the derivatives of interest by the use of (|19bll9c|l . For 
Px ) Qx , Py i Q y ,P u and Q u , we apply d x , d y and d u , re- 
spectively, to <|19bll9c[) and substitute all appearances of 
R u in terms of the other derivatives of interest by the use 
of H19a|) . There is some freedom in choosing the ordering 
of partial derivatives that converts the resulting second- 
order system into a first-order one. The following results 
from one of such choices: 



2d u R x 


— $x-Rx ~r~ $yQx 


+ d x R z , 


(30a) 


2f^u Ry 


— $x-Ry H~ $yQy 


+ d y R z , 


(30b) 


2d u R z 


= d x R x + d y Ry 


+ d z R z , 


(30c) 


d z P x 


= d x R X: 




(30d) 


dzQx 


Oy R X , 




(30e) 




dx Ry j 




(30f) 


d z Q y 


= dyRy, 




(30g) 


2d z P u 


— dx-Rx ~r~ $xQy 


+ d x R z , 


(30h) 


2d z Q u 


— OyP X ~\~ OyQy 


+ d y R z , 


(30i) 



This is again a characteristic problem in canonical form. 
Notice that the first seven equations, l|30a|30g| ) , decouple 
from the last two, since they involve all the variables 
except P u or Q u . Thus R x , R y , R z , P x , P y ,Q x and Q y 
can be found without knowledge of P u or Q u . Once they 
are known, they can be used as known sources in the 
right hand side of I|30hl30i|l to integrate P u and Q u via 

Pu = Pu\z=q + \d x f (P x + Q y + R z ) dz', (31a) 
* Jo 

1 f z 

Qu = Q«|«=o + rfy / (P x + Q y + R z )dz'. (31b) 
1 Jo 

For this reason, we now restrict attention to the subsys- 
tem H30a|30g| ). It constitutes a characteristic problem of 
the form 



B a d a w = 0, 



(32) 



for w = (R x , R v , R Z ,P X , Qx, Py, Qy) where the four prin- 
cipal seven-dimensional matrices B a are symmetric and 
have constant coefficients, which implies a conservation 
law: 



d a (wB a w)=Q. 



(33) 



Integrating this conservation law on the "hyper-prism" 
we have 



w(B % 

u+z=T 



-B z )w dzdxdy 



wB u w dzdxdy 

u=0 



wB z w dudxdy (34) 

z=0 



Since wB u w 
Pi + Ql 



p2 

y 



2{R 2 X + Rl + R 2 Z ) and wB z w = -R\ 
- Q 2 this equation reads explicitly 



f {2Rl + 2Rl + Rl+P 2 x +Ql + Pt + Ql)d-Z T 
= 2 f (Rl + Rl+R 2 z )d^ u 

{~R 2 Z + P 2 x +Ql + P 2 + Q 2 y )d^ z , (35) 
which trivially implies 



w\\% < 2( / (Ri + Ri]+R z z )dX u 



(P 2 + Q 2 x + P 2 + Q 2 v )di: z ), (36) 



Here the norm of the derivatives w is the natural exten- 
sion of the norm of the solution v to seven dimensions, 
defined by 



HI 



(Rl + R 2 V + R 2 Z + Pi + Ql + Pi 



fy)d^T- 

(37) 

The inequality (l.'UiH gives us an estimate of the size of the 
derivatives of interest in terms of the derivatives of the 
free data, accounting for all of the free data without over- 
counting, in a manner similar to the estimates obtained 
in the case of the standard initial value problem, except 
that P u and Q u are treated in a special manner, that is: 
as solutions of ordinary differential equations, Eqs. Ij30hl 
I30i|) . Therefore, the derivatives of the solution of the 
characteristic problem of the wave equation are a priori 
under control. 



IV. WELL-POSEDNESS OF THE 
CHARACTERISTIC PROBLEM OF THE WAVE 
EQUATION WITH A NULL CONE DATA 
SURFACE 

The case of data prescribed on a null cone shares all the 
conceptual features of the case of a null plane. However, 
it becomes technically much more involved, essentially 
due to the presence of non-constant coefficients in the 
expression of the wave equation in spherical coordinates. 
The estimates of the solution in terms of the data can be 
derived in a manner similar as the case of a null plane, 
and are obtained in Subsection llV Al In order to estimate 
the derivatives, however, new tools are necessary, which 
are developed in Subsection IIV Bl 



A. Estimate of the solution 

If we wish to prescribe data on a null cone, it is con- 
venient to transform to spherical coordinates (x, y, z) — > 
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(r, 9, (f>), in which case the wave equation JI} reads 

1 



At ~ -d rr (rip) 1 d g (sin Oip g ) - ., . - 7 

r H sin # r 2 sin 



(38) 



For convenience wc define the new variable 

9 = rip (39) 

and change 6 to a linear coordinate 9 — > s = cos 9, ob- 
taining 



1 



1 



gtt - grr - ^d s ((1 - S 2 ) 9s ) - ———g H = 0. (40) 

Now we change coordinates (t, r, s, 4>) — > (it, r, s, </>) via 

u = t-r, (41) 
and the wave equation (|4()|l becomes 

2g ur -g rr -^d s ((1 - s 2 ). 9s ) - ^ g2) . 900 = 0. (42) 

In order to turn it into a first-order problem we define 
R = g r ,P = g s VT^/r and Q = ^/(vT^sV). We 
obtain the following system 



2i? u — -Rr 



sP 



p, 



n \^l — s 



?Q<P 



r 



Pr 



VT^ U P 

n a 

r r 



1 



R Q 

-tt<h 



g r = R 



(43a) 

(43b) 

(43c) 
(43d) 



For a unique solution to exist in the region of R 4 limited 
by a worldtube of radius r$ and a null cone at u — 0, this 
system requires prescribed values of R on the null surface 
£„ at it = 0, and values of g,P,Q on the worldtube S r 
at r = ro- As in the previous case, the variables P,Q,R 
can be obtained with no knowledge of g, and so can their 
estimate. Ignoring (|43d|) . the remaining equations in the 
system (|43|l have the form 

A u v, u +A r v, r +A s v, s +A+v,i, = Dv (44) 

where v = (R, P, Q), and the matrices are given by 



(45) 



















i 





• 


A r = 1 




1 


• 











V o 









(46) 



A* = - 




, , s/Vj~^l? 
D = — 1 



(47) 



(48) 











1 



There are two difficulties with this characteristic prob- 
lem that are absent in the case of null-plane data of Sec- 
tion ^nj First, the principal matrices have non-constant 
coefficients. This is not an obstacle, though, because, 
since they are symmetric, we can still obtain a "conser- 
vation law". The second difficulty is that there are un- 
differentiated terms, which will appear as sources of the 
conservation law. Multiplying Eq. I|44|l by v on the left, 
and after combining terms appropriately, we have 

(vA n v), u +(vA r v), r +(vA a v), s +(vA*v)^ = 2vDv, 

(49) 

or d a J a = S with J a = vA a v. Here we have 



D 




(50) 



Since vDv = -{P 2 + Q 2 )/r, it follows that 

d a (vA a v) < 0. (51) 

Thus the presence of the undifferentiated terms in this 
particular case will not affect the estimate. We can inte- 
grate now this law on a spacetime volume V , and convert 
the volume integral of the divergence into a surface in- 
tegral over the boundary of V 4 . We choose a region of 
integration limited by E u at it — 0, S r at r = ro and 
at a constant value of it + r = T + ro , with — 1 < s < 1 
and < <fi < 2tt. The region of integration is represented 
in Figure [21 We have 

0> / d a {vA a v)dudrdsd(f> = [ v(A u + A r )v drdsdcf) 
Jv 4 Jet 

vA u v drdsd<j) 



— J vA r v dudsdcf) (52) 

Similarly as in the case of the null plane data surface, we 
define 

|H|| = I (R 2 + P 2 + Q 2 ) drdsd(f> (53) 
= I v(A u + A r )v drdsd<f>, 
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time 




FIG. 2: The region of integration for the characteristic prob- 
lem with a cone data surface. The top boundary surface Et 
lies at u + r = T + ro constant and has measure g(Et = 
drdsd<t>/ y/2. The surface E„ is a null cone at u = and has 
measure dE u = drdsd<f>. The surface E r is a worldtube of 
radius r — ro and has measure dE 2 = dudad(f>. Both surfaces 
E u and E r are data surfaces for the characteristic problem. 

With the inequality reads explicitly 
\\v\\% < 2 J R 2 drdsd<j) 

+ / (-R 2 + P 2 + Q 2 ) dudsd<t> (54) 
which by the same argument as in the case of the null 



plane data surface leads to a similar kind of estimate 




So we see that the "size" of the data on both data surfaces 
controls the "size" of the solution. Finally, once R is 
known, the variable g can be found and estimated by 
integrating (|43d(l with given values on r = r : 

g = g\ ro + f Mr'. (56) 

J To 



B. Estimates of the derivatives 

The results of the previous Subsection establish that 
the size of the solution is controlled by the size of the 
free data. We wish to be able to control the size of the 
derivatives R r , R s , R$, P u , P s , P^, Q u , Q s and Q$ of the 
solution in terms of the derivatives of the free data, as 
well, since these are the derivatives that play a role anal- 
ogous to that of the space derivatives in the case of the 
Cauchy problem. This will be possible if we can write 
down a characteristic system of equations for these nine 
variables, where R r , R s , R^ will be normal variables and 
P u ,P Sl P^, Q Ul Q s and Q$ will be null variables, as a con- 
sequence of the original system of equations. If we write 
down a system for the derivatives themselves in a manner 
completely analogous to the method employed in Sub- 
section |^^BJ the system (|43f) implies the following nine 
equations: 
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2d u R s 



2d u Rd 



-d s P s - 



VT 



r 



1 - s z 
2d u R r ~2 — d s R s 



fd<j,Q s - d s R r 



-y/l — s [ 



--dd>Q<b — drhR r 



1 



r 2 (l — s 2 ) 

d r Ps ~ 

d r P,f, 

d r Qcf, — 



VT^7< 



~VT^7 2 
VT^T 2 



r\/T— s* 



2d r P u - 



VT^s 



2 1 - s 2 1 
d s R r ^—d s P s ^d s 



d^Rr „d<hP s 



ryl — s 



r 2(!_ s 2) 



2sP s sQ -(l + s 2 )P 



sPa, 



r(l — s 2 ) a 



ryl — s 2 



dsRch — d„R r = 



sR s 2V1 - s 2 P s t s{P-rP r )-2Q 4 
r 



d s Rs = -- 



=9,aP s = 



r 

Qs 

r 
p 

dsRch = 

r 

T 



1— 

sRs 



l VT^ 



sP rf 



r(l — s 2 )^ 



rdrkRrh — 



2sP s ( s Q -(l + s 2 )P 2P a 

J 

r 



sP 



2 (l-s 2 ) 



r 



(57a) 
(57b) 
(57c) 
(57d) 
(57e) 

(57f) 
(57g) 
(57h) 
(57i) 



The first seven equations do not involve P u nor 
Q u ; therefore, we may consider Eqs. (|57aJl -( |5Tg| l as 
a system of seven equations for the seven variables 
Rr,R s , R(f>, Ps,P<j>, Qs and the solution of which can 
be used as a known source for l|57hl) - (|57i|) . For the mo- 
ment, we ignore the last two equations, which we will 
come back to once we have an estimate from the first 
seven. The system l|57a|) -( |5Tg| l has many features that 
make it unsuitable for our purposes. In the first place, 
the undifferentiated function P is involved, but it is not 
one of our nine variables. This can be interpreted as a 
non-homogeneous system of equations for the seven first 
derivatives of interest, where P acts as a known forcing 
source. P's presence would definitely affect any estimate 
that may be implied by the system as it stands. Second, 
there are many undifferentiated terms in these equations, 
making it almost certain that they will play a role in any 
estimates. Thirdly, the principal matrices are not sym- 
metric, although their asymmetry is not severe because 
the vanishing coefficients appear symmetrically. As we 
know, the symmetry of the principal matrices would lead 
directly to a conservation law with a non- vanishing source 
of undifferentiated variables, our main goal. 

The first obstacle is a true impediment as far as we 
can see. The standard way to derive estimates for 
non-homogeneous systems of equations is by means of 
Duhamel's principle, which allows one to express the so- 
lution of the non-homogeneous system in terms of the 
solution of the associated homogeneous system and forc- 
ing source function p|. We are not aware of any anal- 
ogous principle for the characteristic problem. Thus, at 
the moment we are forced to consider only homogeneous 



characteristic problems. Therefore, we must find a choice 
of fundamental variables for this system which removes 
the non-homogeneous terms. 

Fortunately, there exists a choice of fundamental vari- 
ables for Eqs. l|57a|) -( |5Tg| l which at the same time sym- 
metrizes the principal matrices and eliminates the ap- 
pearance of the undifferentiated function P, thus taking 
care of the two most important deficiencies of the sys- 
tem of evolution for the derivatives. Instead of using the 
seven first derivatives as fundamental variables, we can 
use the following: 



& = 
Q* = 



Rd 



ryl — s 2 

P/, 



-VT 



r s = 

P S EE 

Q s = 



RsVT 



d s {pVT^) 



QsVT 



(58) 
(59) 
(60) 



This re-scaling will necessarily contribute more undiffer- 
entiated terms to the system because of the dependence 
of the coefficients on r and s. Additionally, the undif- 
ferentiated function P is absorbed into the fundamen- 
tal variable P s . In terms of the variables w = (wi) = 
{R s ,R <t >,R r ,P s ,P' t >,Q s ,Q' 1 '), the system (|57at -(|57g]has 
the form 



B a d a w = Dw, 



(61) 



where all the principal matrices B a are symmetric and 
given explicitly by B" t j = diagonal(2, 2, 2, 0, 0, 0, 0), 
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except 



diagonal(0, 0, — 1, 1, 1, 1, 1) and B s 



tjs r>s rjs 

-D 13 — O 14 — O 26 — 

O 15 — c 23 - « 27 - 



= 

(62) 
(63) 



The matrix D has coefficients that depend on r and s 
but not on the unknown variables. The explicit expres- 
sions of Dij are of no relevance for our purposes because 
they are sufficiently generic to force us to consider the 
most general case. Multiplying by w on the left, Eq. lj6*T|) 
implies 



d a (wB a w) = wDw, 



(64) 



with D — 2D + d a B a . Integrating I^U on the volume 
V enclosed by £y, £„ and S r we find 



jw(B u + B r )wdY, T - fwB u wdY, u - hvB r wdE 



wt)w dV 4 



V 4 



We define 



INI 



d^n 



(65) 



(66) 



with w 2 = J^i w i- Then, exactly as in Subsection IIII Bl 
and because B u + B r is the identity matrix, Eq. 165|) 
implies 

IMIt < 2( jR 2 r +R s2 + R* 2 dZ u 



wDw dV 4 



(67) 



The presence of the undifferentiated terms is a compli- 
cation to the estimate of w in terms of the free data on 
both surfaces. This complication, however, can be re- 
solved in a manner analogous as the case of the initial 
value problem dealt with in Section^J First notice that 



wDw < cw 2 



(68) 



where c = max(|-Dy-|) in the volume V 4 . Thus 
J^wDwdV 4 < c Jq dt f s w 2 dS t . With this, the in- 
equality (|67|l implies 

Nil < 2( fe + i? s2 + i? 02 d£„ 

+ [p s2 + P' t)2 + Q s2 + Q* 2 dT, r 



-c I \\w\\idt 
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(69) 



This holds for any fixed value of T, where E u and S r 
extend as far as their intersection with the surface at 
fixed value of u + r = T. If we denote Yi u ,t and S rit 
the regions of E u and S r extending only so far as their 
intersection with u + r = t < T, we can write 

IHI? < 2 ( [ R*+R s2 +R* 2 di: Uit 



+c / \\w\\ 2 dt' 
Jo 

< 2( /R 2 +P s2 +^ 2 rfS„ 
+ [p^+P^+Q^+Q^dEr 
+c f \\w\\l,dt', 



(70) 



which holds for any value of t < T. Using this inequality 
recursively into the right-hand side of (|69|) we have 



\w\\t < 2(1 + cT 



c 2 T 2 



2 "' ml 



Ips2 + p<i>2 + qs2 + Q<p2 dTir 



+c m+1 [ dh f dt 2 ... [ \Jw\\? dt m+ {71) 
Jo Jo Jo 

for any given non-negative integer m. In the limit for 
ttl — ► oo the sequence in the right-hand side converges if 
cT < 1, in which case we have 



I Mir < 2eCT ( jK + R^ + R^dT.u 

.p9* + Q>* + Q*2,l2 r ) (72) 



This is the final estimate for the derivatives in terms of 
the derivatives of the free data on both surfaces. Like in 
the case of the initial value problem dealt with in Sec- 
tion[n] the estimate involves an exponential factor, essen- 
tially due to the presence of undifferentiated terms in the 
system of equations for w. As usual in such cases, even 
though the estimate is useful in order to prove that the 
solution depends continuously on the data for any value 
of T, it is impractical for large T for the purpose of esti- 
mating the error in a numerical solution. In particular, 
our proof only guarantees the estimate for T < c _1 . Per- 
haps with greater care, possibly by using Gronwall's in- 
equality in integral form, the estimate could be extended 
to longer values of T. 

Since the seven derivatives of interest can be found and 
estimated independently of P u and Q u , they can now be 
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used as known sources for Eqs. i|57h[) - (|57i|l written in 
terms of the new fundamental variables. The solutions 
P u and Q u can be obtained by quadratures as 

Pu = -Pu\r + ^-^d s [ r R r + P s +Q*dr', (73) 

T 1 f T 

Qu = —Qu\ ro + - 7= M Rr+P+Q^dr' '(74) 

r zrvl - s Jr 

from the known functions and from free data given on 
r = ro. 

V. CONCLUDING REMARKS 

Summarizing, Sections IIIII and IIVI develop the proofs 
of the following two theorems: 

Theorem 1 Consider Eqs. \19a\) - \19djl . representing the 
first-order reduction of the wave equation in three spa- 
tial cartesian coordinates (x, y, z) and one null coordinate 
u = t — z. Given data R\ u =q = f(x, y, z) on the null sur- 
face £„ at u — withO < z <T, P\ z= q = g(u,x,y) and 
Q\ z =o — h{u,x,y) on the timelike surface S 2 at z = 
with < u < T, the unique solution v = (i?, P, Q) peri- 
odic in (x, y) satisfies the estimate 

I M It < 2 (Jtfd^u + J(P 2 + Q 2 ) dS 2 ) ■ 

with \\v\\% = J St (R 2 + P 2 + Q 2 )dT, T , where T, T is the 
spacelike surface u + z = T for < u < T and < z < T . 

The derivatives of v are similarly bounded by the 
derivatives off,g and h. 

Theorem 2 Consider Eqs. j^ffo^ - f^cl ), representing the 
first-order reduction of the wave equation in three spatial 
spherical coordinates (r, s — cos 9, 4>) and one null co- 
ordinate u = t — r. Given data R\ u =o — /(r, s,0) on 
the null surface S u at u — with ro < r < T + ro, 
P\ r=ra — g(u,s,<j)) and Q\ r=ro — h(u,s,4>) on the time- 
like surface E r at r — ro with < u < T , the unique 
solution v — (_R, P, Q) satisfies the estimate 

I M It < 2 {jR 2 d^ u + J(P 2 + Q 2 ) dZr^J ■ 

with \\v\\% = J St (R 2 + P 2 + Q 2 )dT, T , where T, T is the 
spacelike surface u + r = T + ro for < u < T and 
r$ < r < T + r . 

The derivatives of v are similarly bounded by the 
derivatives of f,g and h for small values ofT. 

These results are relevant to the stability of the solu- 
tions of the wave equation constructed from data given 
on two intersecting transverse surfaces, one of which is 
timelike and the other one is characteristic. The theo- 
rems guarantee that the solutions will be stable under 



small perturbations of the data on such two surfaces. 
Since the existence and uniqueness of solutions is already 
guaranteed by Duff's theorem 0, our result generalizes 
the standard notion of well-posedness, available in the 
context of Cauchy problems, to the characteristic prob- 
lem of the wave equation. 

In the wider context of general interest in character- 
istic problems of any kind, three publications posterior 
to Duff's pioneering work stand out for relevance 
and motivation. In the first place, Miiller zum Hagen & 
Seifert [6| recognized the value of energy estimates and 
correctly characterize the data surfaces and their role in 
the estimates for all types of problems, including prob- 
lems with one or more characteristic surfaces. Perhaps 
due to the generic nature of their work, their estimates 
appear to have treated indifferently the normal data and 
the null data, making no distinction between free data 
and data that propagates within each characteristic sur- 
faces. This problem is pointed out by Rendall |7|, who, in 
reference to the Miiller zum Hagen & Seifert work, says 
that "they attempted, not entirely successfully, to give 
an existence and uniqueness proof by following step by 
step the treatment of the Cauchy problem." Rendall pro- 
ceeds to detail arguments for existence, uniqueness and 
stability for symmetric hyperbolic systems -like the wave 
equation- in the case of two characteristic transverse data 
surfaces. Rendall does conclude that estimates exist on 
the basis of the well posedness of the associated Cauchy 
problem, but unfortunately, offers no explicit estimates 
of the solution in terms of the set of free data, which 
he identifies. Balean 00,^3 later calculates estimates 
of the energy kind for the wave equation with one char- 
acteristic and one timelike data surface where the null 
data are treated as a flow of information across a time- 
like boundary, and the distinction between the contribu- 
tions of the free data and the transported data is not 
markedly emphasized. All three works deal exclusively 
with second-order equations, leaving Duff's systematic 
approach to first-order characteristic problems without 
a follow-up -in a strictly formal sense. In fact, we have 
been unable to identify any other published literature 
dealing with estimates for characteristic problems of any 
kind. Sometimes the second-order formulation of a prob- 
lem and its associated first-order reduction are regarded 
as equally valid and interchangeable, but we find that 
new insights are often to be gained when a problem (even 
a well-understood one) is viewed from a different vantage 
point, a motivation that underlies our current work. 

The method used here to address the question of sta- 
bility is quite sensitive to the presence of undifferentiated 
terms and is likely to be sensitive to the presence of non- 
linear terms as well. Still, its strength lies in its concep- 
tual features, which depart quite significantly from the 
three predecessors already referred to. It is to be hoped 
that our conceptual framework will be useful as a gen- 
eral guideline for other characteristic problems. Work is 
in progress [TH to extend the method to generic linear 
characteristic problems for first-order systems of equa- 
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